Collective Rydberg excitations of an atomic gas confined in a ring lattice 



00 

o 
o 

(N 
o 

Q 

(N 



> 

00 

oo 

o 



X 



B. Olmos/'B R- Gonzalez-Ferez/'[l| and I. Lesanovsky^'d 

' Instituto 'Carlos F de Fisica Teorica y Computacional and Departamento de Fisica Atomica, 
Molecular y Nuclear, Universidad de Granada, E-18071 Granada, Spain 
^Institute of Theoretical Physics, University of Innsbruck and Institute for Quantum 
Optics and Quantum Information of the Austrian Academy of Sciences, Innsbruck, Austria 

We study the excitation dynamics of Rydberg atoms in a one-dimensional lattice with periodic 
boundary conditions where the atomic Rydberg states are resonantly excited from the electronic 
ground state. Our description of the corresponding dynamics is numerically exact within the perfect 
blockade regime, i.e. no two atoms in a given range can be excited. The time-evolution of the 
mean Rydberg density, density-density correlations as well as entanglement properties are analyzed 
in detail. We demonstrate that the short time dynamics is universal and dominated by quantum 
phenomena, while for larger time the characteristics of the lattice become important and the classical 
features determine the dynamics. The results of the perfect blockade approach are compared to the 
predictions of an effective Hamiltonian which includes the interaction of two neighboring Rydberg 
atoms up to second order perturbation theory. 

PACS numbers: 42.50.Fx,32.80.Qk,32.80.Rm 



I. INTRODUCTION 

During recent years, a new class of experiments in ul- 
tracold gases emerged, dedicated to the study of atoms 
excited to Rydberg states (ij which interact strongly via 
dipolc-dipole or van-der-Waals forces. The most intrigu- 
ing manifestation of this interaction is the Rydberg block- 
ade mechanism Q which prevents the excitation of 
a Rydberg atom in the vicinity of an already excited 
one. On the theory side, this blockade has been thor- 
oughly studied in the context of quantum information 
processing 0, Q for it is a natural implementation of 
a state dependent interaction which is essential to de- 
vise two qubit gates. Moreover, it has been theoretically 
shown that the long-ranged character of the interaction 
can be employed to manipulate whole atomic ensembles 
by just a single control atom Very recently, the state- 
dependent dynamics between two Rydberg atoms, spa- 
tially separated by several micrometers, was observed ex- 
perimentally d, 

In the context of gases, a first experimental indication 
of the strong Rydberg- Rydberg interaction was the non- 
linear behavior of the number of excited atoms as a func- 
tion of increasing laser power and atomic density 0, @] . 
Later, it was shown that, in the case of a dense gas, the 
Rydberg blockade gives rise to the formation of coher- 
ent collective excitations - so-called 'superatoms' [ol. [lo|. 
Due to the strong interaction, the timescale of the excita- 
tion dynamics of a Rydberg gas is typically much shorter 
than those of the external dynamics. As a consequence, 
such gases can be considered 'frozen' [ll|, [l2| in a given 
configuration and the evolution of Rydberg excitations is 
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usually described by a spin model [Tj, [iJI , where the spin 
up/down state represents a Rydberg/ground state atom. 
Unlike in a typical solid state system, there are no signifi- 
cant dissipative processes which make the system assume 
its ground state over the typical experimental timescale. 
Therefore, one can regard the dynamics as fully coher- 
ent [l^ and the time-evolution of quantities like the 
mean number of Rydberg excitations is expected to de- 
pend crucially on the initial state. An experiment, how- 
ever, consists of many successive measurements which 
are performed for different configurations (initial states) 
of the atomic gas, i.e. changing positions of the atoms. 
Hence, the results of many experimental realizations have 
to be averaged in order find expectation values of physi- 
cal quantities. 

In the present work we study the excitation properties 
of a Rydberg gas in a particularly structured and sym- 
metric scenario. In our setup, a large number A^o S> 1 
of ground state atoms per site are homogeneously dis- 
tributed over a ring lattice [17]. At most a single Ry- 
dberg atom per site can be excited via a resonant laser 
that is switched on instantaneously. Due to the underly- 
ing lattice the atomic configuration remains unchanged 
for any experimental realization and no averaging in the 
above sense is required. The atomic motion is consid- 
ered frozen on relevant time scale. In the framework of 
the perfect blockade regime, we study the temporal evo- 
lution of the Rydberg excitation number, the formation 
of correlations in the Rydberg density and the entan- 
glement properties in lattices with up to 25 sites. We 
demonstrate that the dynamics of this system is divided 
into short and long time domains which are, respectively, 
independent and dependent on the lattice size. Our cal- 
culations are numerically exact in the limit of the perfect 
Rydberg blockade and therefore might serve as reference 
for numerical methods that are developed for treating ef- 
fectively one-dimensional many-particle systems. Effects 
that go beyond this regime are treated by second order 
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perturbation theory, i.e. by a Hamiltonian which is ob- 
tained by adiabatically ehminating highly excited energy 
levels. 

The paper is structured as follows: In Sec. [IT] we de- 
rive an effective Hamiltonian which describes the Ryd- 
berg excitation dynamics of atoms confined to a ring lat- 
tice. Section IIIII is dedicated to the discussion of the 
symmetry properties of the system and the consequent 
arising simplification of the numerical calculation. The 
regime of the perfect blockade is thoroughly studied in 
Sec. IIVI In Sec. |V] we discuss the adiabatic elimination 
procedure which we use to account for effects that go be- 
yond the perfect blockade and perform a comparison of 
the results obtained by both approaches. The conclusion 
and outlook are provided in Sec. IVII 



II. THE HAMILTONIAN 

We consider a gas of bosonic atoms that is confined 
to a large spacing optical ring lattice with N sites with 
periodicity a ~ fim. Such lattices have been proposed to 
be created by means of two interfering laser beams with 
A'' being typically of the order of 20 Moreover, there 
is also a way to create a ring lattice approximatively in 
a standard rectangular large spacing optical lattice [l^ 
which is generated by crossed laser beams. To this end 
one can remove atoms from unwanted sites by employ- 
ing the electron beam technique presented in Ref. [l9|, 
thereby 'cutting out' a ring. 

For the sake of simplicity, we assume an uniform and 
large atomic density, i.e. the same number of ground 
state atoms A'o ^ 1 per site. The atoms populate the 
ground state of each lattice site, which for the k-th one is 
described by the Wannier function ^^(x), where x repre- 
sents the spatial coordinates. We assume no hopping and 
hence no particle exchange between the sites. The atoms 
are modeled as a two-level system: the ground state \g) 
and the Rydberg state |r). Experimentally these two lev- 
els are usually coupled by a two-photon transition. Here, 
we assume, without any loss of generality, that they are 
coupled resonantly by a laser of Rabi frequency Qq. The 
corresponding Hamiltonian reads 



N 



fc=i 



(1) 



where bj^ and r|. {bk and r^) represent the creation (an- 
nihilation) of a ground and a Rydberg state atom at the 
fc-th site, respectively. For the sake of simplicity we as- 
sume that the Rydberg atoms experience the same trap- 
ping potential than the ground state atoms. This require- 
ment is, however, not crucial as the typical timescale of 
the electronic excitation dynamics which we are going to 
study is much smaller than the dephasing time due to 
different trapping potentials experienced by \g) and |r). 
Since we are working in the limit of a large number of 
ground state atoms, A'o 3> 1, contained in each site, they 



can be treated as a classical field and we can replace &| 
and bk by \/]Vo. This yields the Hamiltonian 
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{rk+rl), 
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where we assume a relative phase equal to between 
the condensates confined to different sites. Moreover, 
we neglect radiative decay. This is justified because the 
typical timescale of the excitation dynamics is given by 
the inverse of the collective Rabi frequency = fioA/^ 
which can easily exceed tens of MHz. We will later see 
that due to the Rydberg blockade at most N/2 Rydberg 
atoms will be excited on the ring. As a consequence, the 
time for the first emission of a photon scales in the worst 
case like icmit — {'2/N)To, with the single Rydberg atom 
decay rate Tq"^ ~ 100 kHz. Thus, for the ring sizes under 
consideration the condition Q 3> t^^it met. 

Let us now discuss the interaction between the Ry- 
dberg atoms. We focus here on the van-der-Waals in- 
teraction which is given as V{x) = Cq/x^ with x being 
the interparticle distance. It has been shown [lo| that 
even this short-ranged, i.e. quickly decaying, interaction 
can significantly affect the excitation dynamics of atoms 
which are several fim apart. This is rooted in the large 
value of the coefficient Cg , which grows proportional to 
n}^ with n being the principal quantum number of the 



excited level ^20i] . 
by 



The interaction Hamiltonian is given 



^^int = ^ Vkknkirik - 1) + ^y^Vjj 



(3) 



where = rj^rk is the Rydberg particle number operator 
and Vij = / dxdx' |*,(x)|^ |*j(x')|^ V{\x - x'|). We as- 
sume that the spatial extension of the Wannier functions 
^j(x) is much smaller than the lattice spacing a. In this 
case we can write Vij « Ce/{\i — where we make 

the assumption that the range of interactions is much 
smaller than the radius of the ring lattice. Moreover, we 
assume that the on-site interaction is much larger than 
any other in the system, i.e. \Vkk\ ^ This im- 

plies that a double occupancy of a single site is ruled out 
and hence Uk has the two eigenvalues and 1. Within 
these approximations and after introducing the energy 
scale e = hil the Hamiltonian reads 



H — Ho + Hi, 



N 
k=l 



{rk + rl) + ^ AiUkUk+i 



1=1 



(4) 



with A; = Co/[{l a)^ e]. The maximal range of the inter- 
actions is given by m a, with m a positive integer. The 
value of m depends on the features of the considered sys- 
tem, such as the lattice spacing and the strength of the 
interaction. If the interaction rapidly decays as the dis- 
tance is enhanced only a few neighboring sites interact. 
Since for the van-der-Waals potential |V/cfc±2| is 64 times 
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smaller than |Vfefc±i|, our study is focused on the m = 2 
case. 

Without the possibility of a double occupancy, the op- 
erator rj. can be interpreted to create a superatom on the 
A;-th site, i.e. a symmetric superposition of all possible 
single atom excitations on that site. This is done with a 
rate fl. In this spirit, the Hamiltonian ^ describes the 
local dynamical creation and annihilation of such super- 
atoms and their interaction. This Hamiltonian can be 
equivalently formulated as a spin model as it is done in 
Refs. 
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III. SYMMETRIES AND STATE 
REPRESENTATION 

Our goal is to study the dynamics of this system with 
the vacuum state |0) (such that |0) — 0) serving as 
initial state. To that purpose, we have to solve the time- 
dependent Schrodinger equation of the Hamiltonian ^ 
and we perform this task by making use of a basis com- 
posed by all possible configurations in the ring. For an 
increasing site number TV, the exact solution of this prob- 
lem becomes quickly intractable as the dimension of total 
the Hilbert grows as 2^. By exploiting the symmetries 
of the Hamiltonian ^ we can show, however, that the 
dimension of the subspace in which the evolution takes 
place is significantly reduced. 

There are two symmetry operations on a ring lattice 
which are of interest in our system: cyclic shifts by I 
sites and the reversal of the order of the lattice sites. 
The former is represented by the operators Xi with I — 
1,2,...,N where Xi = X[, while the latter operation 
is the parity and is denoted by R. The action of these 
operations on the creation and annihilation operators is 
defined through 

X-\kXi=rk+i X-\iXi=ri^, (5) 

R-^rkR^rN-k+i R'^Ar ^ r^^_^^^, (6) 

from which follows that R and all Xi are unitary, i.e. 

= i?t and X^^ = Xj. The Hamiltonian Q is in- 
variant under these operations, i.e. [H,Xi] = [H,R] — 0. 

Let us now focus on the vacuum state |0), where no 
Rydberg atom is excited. This is an eigenstate of all 
cyclic shifts and the reversal operator with eigenvalue 1: 

Xi\0) = \0} i?|0>=.|0). (7) 

Only a small subset of the 2^ states, spanning the whole 
Hilbert space, actually has these properties. Each of 
these maximally symmetric states can be understood as 
a superposition of all states that are equivalent under ro- 
tation and reversal of the sites. Since the Hamiltonian 
conserves the symmetries of the initial state, the evo- 
lution of the system will take place in the subspace of 
the Hilbert space spanned by the maximally symmetric 
states. By using these states, the dimension of the prob- 
lem is dramatically reduced. For example, for N = 10 the 



FIG. 1: Energy level structure of the Hamiltonian Q for 
m — 2. The spectrum consists of highly degenerate subspaces 
which are labeled by i/ energetically separated by eA. The 
laser {Ho) causes an energy splitting of the degenerate levels. 
In addition, it couples states belonging to a given i^-subspace 
and connects subspaces with \v — v'\ — 1, 2. 

dimension of the basis decreases from 2^" — 1024 down 
to 78. For our computations we need an algorithm to 
quickly generate the maximally symmetric states among 
which the evolution takes place. Such an algorithm is 
presented in Ref. [23|. There, these states are called 
bracelets, and are recursively generated in an optimal 
way. The amount of CPU time grows only proportional 
to the number of bracelets produced. 



IV. PERFECT BLOCKADE 

We will consider from now on the case m = 2 if not 
explicitly said otherwise, i.e. only neighboring sites in- 
teract. The Hamiltonian ^ is then governed by two 
energy scales: the one associated to the laser excitation 
Hq, i.e. the collective Rabi frequency e — h^} and the one 
related to the Rydberg- Rydberg interaction Hi^t, given 
by eA = eAi. 

The spectrum of the interaction Hamiltonian iJjnt de- 
composes into + 1 degenerate subspaces of energy 
— v eA with ly 0,...,N counting the number of 
pairs of neighboring excitations. The laser Hamiltonian 
Hq is switched on instantaneously and drives the dynam- 
ics within a given j/-subspace and couples subspaces with 
\i/ — < 2 (see Fig. [T]). The timescale associated with 
the evolution inside a i^-subspace is tq = h/e, whereas 
the typical time of inter-subspace transitions is given by 
Tint — fiA/e = To A. We consider here the regime in 
which A ^ 1, i.e., the interaction energy of two neigh- 
boring Rydberg atoms is much larger than the collective 
Rabi frequency and hence tq ^ Tint- For the van-der- 
Waals case, this condition is hil <C jCgl/a^. 

The physical initial state is the vacuum, |0) which be- 
longs to the subspace with v — 0. In our approxima- 
tion we neglect the coupling between this subspace and 
those including higher excitations. This is the perfect 
blockade approach, which is valid for t <^ Tint, i-e. the 
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O-exc 1-exc 2-exc 3-exc 4-exc 5-exc 



FIG. 2: Graph with the basis of states for = 10 in which 
the time-evolution takes place. In each column a subspace 
of a given number of Rydberg excitations is shown, see text 
for further information. The laser (Hamiltonian Ho) couples 
only states belonging to adjacent subspaces. The coupling 
strength (transition probability) between the individual states 
is encoded in the colors. 



time it takes to perform a transition between adjacent 
i/-subspaces. The restriction to the 1^ = subspace leads 
to a further reduction of the dimension of Hilbert space 
in which the temporal evolution takes place. For exam- 
ple, for = 10 the number of states to be considered in 
the basis set expansion decreases from 78 to 14. This is 
to be compared to the 1024 states which span the entire 
Hilbert space of the system. 

The laser Hamiltonian Hq couples states in the Hilbert 
space whose number of excitations differ by one. The 
corresponding maximally symmetric states and the cou- 
pling between them by means of Hq can be graphically 
illustrated as shown in Fig. [51 The way these states are 
coupled is qualitatively similar for different lattice sizes 
so, for simplicity reasons, we discuss here the lattice with 
N = 10. In Fig. [2] the states are denoted by the number 
of Rydberg excitation, and a subscript is added when 
more than one configuration with the same excitation 
number is possible. Note that, for an even (odd) number 
of sites the maximal number of Rydberg atoms in one of 
these states is ^ {^^y^), e.g. |5) for 10 sites. 

Starting from the vacuum, there are several excitations 
paths with different probabilities that connect the states. 
The larger the amount of Rydberg atoms, the more con- 
straints are found to allocate the next excitation. As a 
consequence, we encounter several excitation paths that 
do not reach the state with the maximal number of Ry- 
dberg excitations, but end in others, such as the |4b) 
and |4c) states for N = 10. The features of these frus- 
trated states strongly depend on the lattice size, and their 
amount increases as N is increased. In particular, their 
existence provokes quantitative differences in the dynam- 
ics of two lattices with different A^- value for large times. 



This is reflected in the time-evolution of all quantities we 
are going to study throughout this work. The dynamics 
can be always divided into two different domains. For 
short times, t < 4to the behavior is universal, i.e., in- 
dependent of the size of the lattice, whereas for longer 
times a dependence on N is observed. 

A. Two-sites density matrix 

The reduced density matrix of two neighboring sites 
is needed to investigate the temporal evolution of local 
properties such as the mean density of Rydberg atoms or 
the entanglement between two adjacent sites. Since the 
wavefunction |^'(t)) is spanned in the subspace of fully 
symmetric states, all sites are indistinguishable and we 
can take 1 and 2 as representative adjacent lattice sites. 
The two-sites reduced density matrix is obtained from 
the density matrix of the full system p{t) = \^ (t)) {"^i {t)\ 
by performing the partial trace over all the remaining 
sites, 

p(i2)(i)^Tr3.4 ^(p(i)). (8) 

The basis for the two-sites states is {\gg), \gr), \rg), \rr)}. 
The restriction to the fully symmetric subspace imposes 
{gr\p<^'^^\gr) = {rg\p^'''>\rg) ^ P and {gg\p<^'^^gr) = 
(f/g|p^^^^ |rg) = 7, while the perfect blockade prevents 
the excitation of atoms in two neighboring sites and 
hence the entries {A\p^^^^\rr) and (rrlp'^^^^A) are zero 
for \A) = {\gg),\gr),\rg)}. As a consequence, the re- 
duced density matrix has the particularly simple form: 
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where a, (3, 7 and S are four time-dependent (complex) 
parameters. In particular, a and f3 are real, and due to 
the normalization of the wavefunction, Tri^2 P*^^^-* (0 — 
a + 2f] = 1. Hence, only three of these parameters are 
independent. Performing the trace over the states of site 
number 2 yields the single particle density matrix 

p(^)(0=Tr.(p(-)(t)) = (l-^ (10) 

B. Rydberg density 

The first local property under consideration is the time 
evolution of the expectation value of the Rydberg density 
rifc. By using Eq. PH)) one obtains 

(nfc(t)) -Tr (t)nfe(t)) = /3, (11) 

and the total number of Rydberg atoms evaluates to 
(iVRyd(i)) = N (ukit)) = Np. In Fig. Hi we show {ukit)) 
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FIG. 3: Expectation value of the Rydberg density versus 
time for (a) four different ring sizes and t < 100 ro, and (b) 
detail of the short time evolution for A*' = 10 and 25. The 
computations have been performed assuming perfect blockade 
of the adjacent neighbor. 



FIG. 4: Short time behavior of the density-density correla- 
tions g2{k,t) for A'^ = 25. Correlations emerge successively 
during the time evolution. Due to the perfect blockade con- 
dition strong oscillations of g2{k,t) with a period k — 2 are 
observed. 



C. Density-density correlation function 



as a function of time for different lattice sizes. In the 
right panel (Fig. [Sb) a magnified view of the short time 
dynamics is provided for = 10 and 25. 

We observe at first a steep increase, which is propor- 
tional to , that culminates in a pronounced peak located 
a.t t = 1.09 tq. This peak is independent of the lattice 
size, as we are still witnessing the short time behavior. 
For much larger times, (nfc(i)) becomes dependent of N 
and oscillates with a frequency / ~ 0.48 fl about a mean 
value of {nk{t)) « 0.26. This mean value and also / turn 
out to be independent of the ring size, however, the exact 
shape of {n]^{t)) strongly depends on N . A similar result 
and a possible explanation of the nature of this value 
is given in Ref. |13j . The time averages are performed 
by means of numerical integration over time in a large 
enough interval [5ro,200ro]. The lower limit of this in- 
terval is chosen large enough to avoid the initial effects of 
turning on the laser. The higher one is kept shorter than 
the corresponding revival time. The amplitude of the os- 
cillations decreases considerably with increasing lattice 
size. This amplitude can be measured by means of the 
standard deviation. For iV = 10, a quasi-steady state 
with large fluctuations characterized by a standard de- 
viation u{{nk(t))) = 0.062 about {nk{t)) is established, 
while in the case of TV = 25 these fluctuations become 
smaller cr((nfc(t))) = 0.0065. 

We have also p erform ed calculations of time-averaged 
Rydberg density {nk{t)) for m = 3 and to = 4 and N > 

0.17 



20. The results for the two cases are (nfc(i))m=3 
and (?T.fc(i))m=4 ~ 0.12, which is an indication for the 
scaling {nk{t))n, ~ (2to)"^ 



The equal-time density- density correlation function 
between two sites denoted by i and j separated by a 
distance |i — j| a = fc a is given by 



(12) 



It measures the conditional probability of finding an ex- 
cited atom at a distance k a from an already excited one 
normalized to the probability of uncorrelated excitation. 

Figure m illustrates the initial evolution of g2{k,t) in 
the time interval [0, 5 tq] ior a, N — 25 lattice. Due to 
the perfect blockade condition, (72(1, i) = for any time. 
The temporal and spatial structure can be understood 
by observing the properties of the laser Hamiltonian Hq 
which drives the excitation dynamics. At the beginning 
only a single particle at site 1 is excited, so the probabil- 
ity of excitation of a second atom is uniform, it can occur 
at arbitrary position (except at a distance of fc = 1). As 
a consequence, there are no correlations for very short 
times, and they emerge successively as time increases. 
The augment occurs at /c = 2, and the density-density 
correlation function for that distance reaches a maximum 
at t « 1.5 To . The high probability of finding two excita- 
tions at the distance fc = 2, i.e. a large value of 52(2, t), 
automatically gives rise to a decrease of 172(3, t) due to 
Rydberg blockade. For larger times, a regular pattern 
of enhanced and suppressed density-density correlations 
characterizes the dynamics. The regular pattern of the 
density-density correlation functions at short times is lost 
as time increases. Here, 172(^7 1) exhibits pronounced fiuc- 
tuations and rapid oscillations around a mean value for 
t >4to. 

In Fig. O we show the time-averaged density-density 
correlation function, (72(^1 1), in the stationary long time 
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FIG. 5: Time-averaged density-density correlations g2(k,t) 
for A'' = 25. For dist ance k — 1, the value is zero due to 
the perfect blockade. g2{k,t) assumes a maximum for k = 
2 (next-nearest neighbor). For larger distances, only weak 
correlations are visible. 
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FIG. 6: (a) Time-evolution of the two-party correlation mea- 
sure Mc (p'^'^^) ^'^^ long times and various lattice sizes, (b) 

Short time behavior of Mc ^p'^^'^ and its classical counter- 
part for iV = 25. 



regime as a function of k, for a lattice of 25 sites. This 
function shows a maximum for fc = 2, while for larger 
intersite distances it approaches the constant value 1, i.e., 
no correlations. As a consequence, we conclude that the 
density-density correlations are only short ranged after 
the initial period in which also long ranged correlations 
are of importance. 

D. Entanglement 

We study the quantum and classical correlations and 
the entanglement of two neighboring sites in this system 
by means of the two-party correlation measure [22l | and 
the entanglement of formation [1^. These quantities can 
be directly related to the entries of the reduced density 
matrix discussed previously. 

Two-party correlation. The two-party correlation 
measure [22] is based on the trace distance |^ and it 
is defined as 

A/c(p(^^))=^Tr|p(i2)-pW®p(2)|^ (13) 

where |^| = VaTA is the positive square root of A^A. Its 
physical meaning is the distance between the state p*-^^^ 
and its reduced product state p^^^ (g) p^^^ It takes into 
account both the classical correlation between two sites 
and the quantum coherence. It generalizes the classical 
distance in the sense that if the two operators commute 
then it is equal to the classical trace or Kolmogorov dis- 
tance between the eigenvalues of p*^^-^^ and p'^^-' ® p'^^K 
We show the time evolution of this correlation measure 
in Fig. [S^ for different sizes of the ring. Initially, for 
the vacuum state, there are no correlations. Analogously 
to the previously analyzed quantities, Mq exhibits an 
A'^-independent short time behavior which here is char- 
acterized by large amplitude oscillations. It is followed 



by an A^-dependent regime, where Mc presents smooth 
oscillations around the mean value Mc — 0.19. As ex- 
pected, the amplitude of these oscillations decreases with 
increasing lattice size. 

We are now interested in finding a classical counterpart 
to this correlation measure. To this end, we make use of 
the density matrix properties. The diagonal elements of 
a density matrix represent the probability of finding the 
corresponding configuration of the sites. For example, in 
p^-^^\ (see equation a and P represent the probabil- 
ity of the two sites to be in the states \gg) and \gr) or 
\rg), respectively; note that the sites are indistinguish- 
able. In the same way, the diagonal components of the 
reduced product density matrix p*^^^ (g) p^^^ provides the 
probability of the two sites being in the corresponding 
product state, e.g. (1 — /3)^ for the state \g) (Xi \g)- 

We take these diagonal elements d[^^^ and d^^^^^ of 
the matrices p^^^-* and p*-^-* (8) p^^\ respectively, as classi- 
cal probability distributions. The Kolmogorov distance 
between these distributions, is defined here as: 

^cla.s(^(12)^^2^|,(12)_^(l«2)|^ (14) 
i—1 

and it provides a classical measure of the two-party corre- 
lation. In terms of the parameters of the density matrix 
this quantity is reduced to 

The classical and the total two-party correlation func- 
tions are presented in Fig. [6)3 for t < 4ro and 25 sites. 
One of the main features due to the quantum behavior 
of the system is the appearance of the two consecutive 
peaks of Mc att = 0.88 tq and i = 1.09 tq. Note that the 
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classical counterpart M^'^^ reproduces only the second 
maximum. Hence, the existence of the first one can be 
only justified by quantum arguments. Due to the abso- 
lute values in expression (jl3p . two discontinuities appear 
in the derivative of Ale around t w 1.9 tq and t w 2.3 tq. 
They are, however, not observed in Mq'^^^ since it only 
depends on the single, smoothly varying, parameter of 
the density matrix, /3. 

To get a deeper insight into the quantum effects on the 
correlations, the difference between the total two-party 
correlation and the classical measure is shown in Fig. [7| 
as a function of time. We have performed a fit to the local 
maxima of this numerical difference using an exponential 
decreasing function. The contribution of the quantum 
correlations loses importance as time is increased and, at 
the same time, the classical dynamics starts to dominate 
the correlations between two neighboring sites. 




FIG. 8: Time-evolution of the entanglement of formation for 
(a) long time and several sizes of the system and (b) short 
time for N = 25. 



— N=25 
---Fitting exp 



F(th0.2378 exp(-0.6767 1) 




FIG. 7: Difference between the two-party correlation mea- 
sure and its classical counterpart for N = 25. The dashed 
line corresponds to an exponential fit to the envelope of this 
difference. 



Concurrence and entanglement of formation. For a 
general state of two qubits represented by means of its 
two-particles density matrix p, the concurrence is given 
by m 



C{p) = max {0, Ai - A2 - A3 - A4} , 



(16) 



where the Ai are the square roots of the eigenvalues, in 
decreasing order, of the matrix pp, where p is the flipped 
matrix of the two-qubit general state p, i.e.. 



{ay ® Uy) (^P^^^^y [Oy ® Gy) . (l7) 



The entanglement of formation of a state of two qubits 
is defined as 



E{p) 



1 + v/1 - C{pf 



(18) 



with /i(x) = — xlog2 a;— (1 — x) log2(l — a;). This quantity 
provides a measure of the resources needed to create a 
certain entangled state, and its range goes from to 1. 



Using the two-sites density matrix describing our sys- 
tem (|9]), and its corresponding flipped matrix, we obtain 
the following A^: 



Ai = (/3+| 



A2 = |/?- 



A3 = A4 = 0. (19) 



These values give rise to two different regimes for the 
concurrence: 



C 



2|<5| 
2/3 



P>\5\ 
P<\5\ 



(20) 



The first condition (3 > \5\ always holds for any size N 
of the lattice, so the concurrence yields 



C 



=2\5\. 



(21) 



The time evolution of the entanglement for the lattices 
with sites N = 15, 20 and 25 is presented in Fig. [8^, and 
an enhancement of the behavior at short times for iV = 25 
is shown in Fig. ^p. Again, two different time domains 
can be distinguished. For short times, the entanglement 
of formation is independent of the ring size. Its maximal 
value, E [p) — 0.23, is reached ai t = 0.73 tq; for a fur- 
ther increase of time, E [p] drastically decreases, e.g. the 
second peak ait — 2.05 tq is reduced roughly by 80%. In 
the long time regime, the entanglement becomes weaker 
with E [p) eventually approaching zero with character- 
istic fluctuations for each N . The amplitudes of these 
fluctuations become smaller as N is increased. 



V. NON-PERFECT BLOCKADE 

The above discussed phenomena have been investi- 
gated in the perfect blockade regime, which assumes that 
the energy scale associated to the van-der-Waals interac- 
tion is infinitely large compared to the one related to the 
laser interaction, i.e. A ^ 1. Since the initial state has 
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no Rydberg excitations, the dynamics of the system is 
restricted to a small subspace including those eigenstates 
of Hint with eigenvalue zero. However, for finite values of 
A, the states with i/ = are coupled to those with ly > 
and, even more, these higher excitations might influence 
the dynamics. 

In this section, we go beyond the perfect blockade ap- 
proach and explore these couplings including their ef- 
fect of up to order 1/A. We thereby derive an effective 
Hamiltonian by dividing the eigenstates of Hint into two 
sets, characterized by their respective quantum number 
ly. The first set of states is formed by the subspace ly — 0, 
whereas the second one contains the rest of energetically 
high- lying excitations. In this framework, the Hamilto- 
nian can be written as 



/ PHP PHQ \ 
^ - ' QHP QHQ J 



(22) 



where P and Q are the projectors on the subspaces with 
energy v = and ly > 0, respectively. A general wave- 
function can be decomposed as 



(23) 



and the time-dependent Schrodinger equation reads: 
^^dt ( ^„ ] = yQHP QHQ ) I n,T, ) ■ (24) 



Due to the large energetic gap between the different sub- 
spaces and since the initial state is the vacuum, the tran- 
sition probability to states with ^ Eq is very small. 
Thus, we can introduce an approximation assuming that 
the time variation of Q"^ is very small and can therefore 
be neglected, i.e. 9t(Q^) = 0. Hence, the equation of 
motion (12411 is reduced to 



ift(P*) = {PHP - PHQ{QHQ)-^QHP) [P^). (25) 

Note that, in this expression, PHP — Hq is the Hamil- 
tonian within the perfect blockade regime. Whereas, the 
second term provides the first correction to this Hamil- 
tonian and represents the contribution of the couplings 
between the v — and v > subspaces. In practice. 
Ho only couples the v = Q subspace and those with 
V — land 2, see Fig. [TJ As a consequence, the Hamilto- 
nian can be rewritten as 



H = 




(26) 



where the subscript R denotes the energetic levels with 
i/ > 2, and the f2ab represent the part of the Hamiltonian 
that couples the states of the subspaces with v = a and 
v' = h. In expression P^ . QHQ can be decomposed into 
the sum of a diagonal matrix. A, including the interaction 



45 p 

40- 

35- 
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25- 
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Energy (e) 



_J \ I L 
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Energy (e) 



FIG. 9: Histogram of all the eigenvalues (density of states) 
of the full Hamiltonian for a system with N = lb, m = 2 
and A = 20. The parameters are chosen such that the indi- 
vidual i/-manifoIds are still recognizable. The inset shows a 
magnified view of the manifolds with v = 2, which are 
broadened by the interaction with the laser. 



between the subspaces, and a full matrix containing the 
couplings fj. 



QHQ 



eA 
2eA 
Afl 



^Rl ^R2 ^R 



A + n. 



The inverse of this matrix can be approximated by 
1 



{QHQ)-' 



A + n 



(27) 



(28) 



where we have used the Neumann series (I — T) = 
X]^o P" ^ square matrix T whose norm satisfies that 
||T|| < I. Since A ^ flab for any a and 6, this condition 
is accomplished for T = A~Tt. Finally, we obtain the 
following expression for the effective Hamiltonian 



Hcff — Ho 
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eA 



2eA 



+ 0(1/A2), (29) 



where we only consider the first three terms and neglect 
higher order corrections. 

Let us now discuss the regime of validity of the ap- 
proximate Hamiltonian (|29p . To this end it is instructive 
to study a case in which the full Hamiltonian ^ is nu- 
merically tractable. This however, can only be done for 
a small number of sites. In the absence of the laser the 
eigenstates of Hamiltonian ([4]) are those of iJint, i.e. the 
highly degenerate i/-manifolds. As soon as the laser is 
turned on this degeneracy is lifted and all the v manifolds 
split up. However, if A is sufficiently large the manifolds 
are still well separated. This regime is presented in Fig. [9] 
where we show a histogram of the eigenenergies (density 
of states) for a lattice with iV = 15 and A = 20. Since 
in this case the system can contain at most 15 pairs of 
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FIG. 10: The Rydberg density as a function of time for = 
20 computed with the perfect blockade treatment, and with 
the adiabatic ehmination scheme using A = 25 and 35. 



FIG. 11: Correlation function for fc = 2 as a function 
of time for = 20 computed with the perfect blockade 
treatment, and with the adiabatic elimination scheme using 
A = 25 and 35. 



consecutive Rydberg atoms, we observe 16 manifolds, i.e. 
< < 15. The energetic separation between the cen- 
tral states of two neighboring subspaces is given by A. 
A magnified view of the spectral structure for the low- 
lying excitations is shown in the inset of Fig. [S] Within 
the framework of the adiabatic elimination the contri- 
bution of the v — 1 and 2 manifolds is included up to 
order ^ in the effective Hamiltonian (|29p . The validity 
of this approximation is restricted to parameter regimes 
in which states belonging to different manifolds are en- 
ergetically well-separated, e.g. two adjacent manifolds 
must not overlap. For = 15, A = 20 is the minimal 
value needed to ensure this separation. For larger lat- 
tices sizes, the value of A has to be increased since with 
growing TV the iz-manifolds contain more and more states 
and thus become successively broader. For example, the 
width of the v = {) manifold scales proportional to N . 

We have investigated the dynamics of a ring with 
iV = 20 sites in the framework of the adiabatic elimi- 
nation using A = 25 and 35. In Figs. [TO] and [11] we 
show the Rydberg density and a density-density corre- 
lation function (for fc = 2) and compare them with the 
results obtained within the perfect blockade approxima- 
tion. The occurring deviations are small. Only minor 
differences are observed at large times, and we encounter 
relative errors below 6.5% and 4% for A = 25 and 35, 
respectively. The results show that the approximated 
inclusion of higher i^-subspaces in the dynamics does only 
lead to small quantitative changes in the behavior of the 
investigated quantities. As anticipated, the deviations 
reduce significantly as A is increased. More qualitative 
differences are expected to occur if the r~^-tail of the 
Rydberg- Rydberg interaction is properly accounted for. 



VI. CONCLUSIONS AND OUTLOOK 

In this work we have performed a numerical analysis 
of the laser-driven Rydberg excitation dynamics of atoms 
confined to a ring lattice. By exploiting the symmetry 



properties of the system and employing the assumption of 
a perfect Rydberg blockade we were able to perform nu- 
merically exact calculations in lattices with up to N=25 
sites. Our findings show that the temporal evolutions of 
the physical quantities, e.g. the Rydberg density and the 
density-density correlations, can be divided into two do- 
mains. For short times, one observes an N-independent 
universal behavior with large amplitude oscillations. For 
longer times, the dynamics is crucially determined by the 
lattice size and the analyzed quantities appear to assume 
a quasi steady state with only small temporal fluctua- 
tions. Moreover, we studied the evolution of the entan- 
glement as well as the quantum and classical correlation 
of two neighboring sites. By separating the quantum 
and classical part of the two-party correlation we showed 
that quantum correlations between neighboring sites de- 
cay rapidly as time passes. In addition, the entanglement 
between neighboring sites turned out to be weak in the 
long time limit after a quick initial increase. 

We eventually relaxed the perfect blockade condition 
by taking into account higher excitation subspaces via 
adiabatic elimination. Propagating the initial vacuum 
state with the corresponding effective Hamiltonian has 
only small effect on the time evolution of the investigated 
quantities. More quantitative changes are expected when 
including the long-ranged tails of the interatomic inter- 
action potential. This however requires more powerful 
numerical methods, such as t-DMRG [1^, which goes 
beyond the scope of this work. 

In the present work we have been focusing on the dy- 
namical properties of this system. A next step would be 
to investigate the corresponding static properties, such 
as eigenstates and eigenvalues. However, since physically 
the system is at i = in the vacuum state it remains 
an open question how the eigenstates can be actually ac- 
cessed, for instance via an adiabatic passage incorporat- 
ing a time-dependent change of the laser detuning and its 
Rabi frequency. In addition, a rather natural extension 
would be to analyze the dynamics of these many-particle 
systems by means of a two dimensional description. Cer- 
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tainly, it is also of interest to explore lattices with dif- 
ferent geometries, e.g. square or triangle. The different 
underlying symmetries arc expected to significantly affect 
the time evolution of these systems. 
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